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Abstract. 

The critical properties of the stochastic susceptible-exposed-infected model on a 
square lattice is studied by numerical simulations and by the use of scaling relations. 
In the presence of an infected individual, a susceptible becomes either infected or 
exposed. Once infected or exposed, the individual remains forever in this state. The 
stationary properties are shown to be the same as those of isotropic percolation so that 
the critical behavior puts the model into the universality class of dynamic percolation. 


PACS numbers: 05.40.-a, 02.50.Ey, 87.10.Hk 


Critical properties of the susceptible-exposed-infected model on a square lattice 2 

1. Introduction 

Spatio-temporal structures as well as fluctuations are essential features of epidemic 
spreading [DSEli]. A description of epidemic spreading that takes into account these 
essential features is provided by stochastic lattice models [siEiEiEiEiiininiiniiisi 
[Til IT5| [T6] in which each site of a lattice is occupied by an individual that can be in 
one of a certain number of states. In the susceptible-infected-recovered (SIR) model 
[I311I111I5], an important model in this context, the possible states are susceptible (S), 
infected (I) or recovered (R). The SIR model is composed by two processes. One in 
which a susceptible becomes infected by a catalytic reaction, S-|-I^I-|-I, and another in 
which an infected becomes recovered spontaneously, I—?-R. Another model, the one that 
will be the object our study here, is the susceptible-exposed-infected (SEI), introduced 
by Tome and de Oliveira [15], in which each individual can be susceptible (S), exposed 
(E) or infected (I). This model has also two processes. In the presence of an infected 
individual, a susceptible individual may become either infected or exposed, processes 
represented by the reactions S-l-I—and S-l-I—)-E-|-I, respectively. 

The distinguish features of these two models are as follows. The spreading of 
the epidemic occurs as long as there are active sites, which are the sites occupied 
by a susceptible individual next to an infected individual. When the active sites 
have disappeared, the infection reaction S-l-I—no longer takes place, that is, no 
new infected individuals are created, and the whole process eventually stops. The 
system Ends itself in one of many absorbing states which, in the SIR model, are the 
conhgurations formed by R and S sites and, in the SEI model, are the conhgurations 
without any pair of neighboring SI sites. Starting from a conhguration full of susceptible 
individuals except for a single infected site, the process generates a cluster of inactive 
sites, which are the R sites of the SIR model and the I sites of the SEI model. In the 
stationary state, which is an absorbing state, these clusters have the same properties as 
the clusters occurring in isotropic percolation model so that the stationary properties 
of these two models are similar to percolation. 

When the rate of infection is small there is no epidemic spreading, that is, the 
generated clusters are all hnite. If however the infection rate is large enough, an 
inhnite cluster of inactive sites is generated and the epidemic spreading takes place. 
The transition from non-spreading to spreading is regarded as a continuous phase 
transition with critical behavior within the universality class of dynamical percolation 
P [ini Hn [121 ESI EH EZl eh ESI [ 20 ]. in fact, the clusters generated by the rules of 
the SEI model can be exactly mapped into the clusters of site percolation so that the 
stationary properties of the SEI model are identical to the properties of site percolation. 
In this sense, it is similar to the model introduced by Alexandrowicz [21] to generate 
percolation clusters. The SEI model can thus be regarded as a standard example 
of a model belonging to the dynamical percolation universality class. Here we are 
concerned with the critical behavior of the SEI model on a square lattice, particularly 
with the numerical calculation of the dynamic critical exponents. As to the static critical 
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exponents, they are the same as the percolation in two dimensions and known exactly 
|22] . Here we perform calculation for the SEI model with results for the dynamic critical 
exponents that are very accurate, and in agreement with the exponents of the dynamic 
percolation universality class [T3l 120] , 

2. Model 

The SEI model is a continuous time stochastic markovian process dehned on a lattice 
where each site can be in one of three states: S, I or E. The allowed transitions are 
those in which just one site changes its state. The transition rate of the process S—j-I is 
bf, where / is the fraction of I sites in the neighborhood of the site to be updated and 
b is the infection parameter. The transition rate of the process S—)-E is af, where a is 
the exposition parameter. Other transitions are forbidden so that sites in states I or E 
remain forever in these states. For convenience, we dehne the parameters p = b/{a -\- b) 
and g = a/(a + 6) so that p -\- q = 1. The model displays two regimes. One in which 
there is no epidemic spreading, occurring for small infection rate, and the other in which 
the epidemic spreading takes place, occurring for large infection rate. The stationary 
properties are close related to site percolation model. In fact, the cluster of I sites can 
be exactly mapped into the cluster of site percolation in the same lattice, in which 
each site is occupied with probability p, as will be shown below. It is well known that 
the percolation model [22] shows a phase transition from a state with hnite clusters, 
occurring for p < p^, to a state with an inhnite cluster, or percolating state, for p > p^, 
where pc is the critical concentration. We thus expect for the SEI model, a transition 
from a non-spreading regime, for p < pc, to a spreading regime, for p > p^. 

The simulation of the model is carried out as follows. Each site of a regular lattice 
with N sites can be in one of three states: occupied by a susceptible (S), by an infected 
(I) or by an exposed (E) individual. At each time step a site is chosen at random. If it 
is in the I or E states nothing happens. If it is in the state S then with probability pf it 
becomes I and with probability qf it becomes E, where / is the fraction of I sites in the 
neighborhood of the chosen site and p + q = 1. Equivalently, if the chosen site is in state 
S, we may randomly choose a neighboring site; if it is in the I state then the chosen site 
becomes I with probability p or E with the complementary probability g = 1 — p. The 
time is increased by an amount equal to 1/N. The neighborhood of a site is dehned as 
its nearest neighbor sites. 

An alternative approach, useful for time dependent simulations, is carried out as 
follows. At each time step, a site is chosen at random from a list comprising the I sites 
only. Next, one of its neighbors is chosen at random. If this neighbor is in state I or E 
nothing happens. If it is in state S, then it changes to I with probability p or to E with 
probability g = 1 — p. The time is increased by an amount equal to l/nj where nj is 
the number of sites in the list, which is the total number of I sites in the lattice, and the 
list is updated. An even more efficient algorithm is set up by using a list of active pair 
of sites, more precisely, a list of nearest neighbor pairs of susceptible-infected sites. At 
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each time step, a pair of the list is chosen at random and the S site of the pair becomes I 
with probability p or E with probability q = 1 — p. The time is increased by an amount 
1/nsi where nsi is the number of entries in the list, which is the total number of SI 
pairs in the lattice, and the list is updated. 

Let us consider a hnite lattice full of susceptible individual except for one site 
which is occupied by an infected individual. The system evolves in time and eventually 
reaches its hnal state which comprises a connected cluster of infected sites in addition 
to exposed and susceptible sites. The exposed sites are found at the boundary of the 
cluster separating the infected sites from the susceptible sites, so that active sites are 
absent. The clusters generated by the SEI rules are the same clusters of site percolation, 
as we show next, in which the sites are occupied independently with probability p. 

The mapping of the stationary properties of the SEI model into the site percolation 
model can be understood as follows. Suppose that a cluster C of infected sites of a 
lattice has been generated by one of the algorithm above, starting from an infected site 
at the origin. Suppose moreover that the dynamics has come to a halt so that there are 
no pairs of the type SI. The boundary B of this cluster is therefore composed only by E 
sites. During the dynamics, whenever the site i becomes either I or E, we keep the used 
random number It is clear that, if site i has turned into an I site then < p; if site i 
has turned into an E site then > p. Let us now consider a replica of the lattice, with 
all sites empty except the site at the origin which is occupied. Next the site i of the 
replica is occupied if fi < p and remains empty if fi > p. By means of this procedure, 
which is the procedure used in site percolation, a cluster of occupied sites is generated, 
which is thus identical to the cluster of I sites of the original lattice. 

3. Scaling relations 

Around the critical point, the quantities that characterize the system are assumed to 
obey a scaling relation. We assume two types of scaling relation. The hrst one is a 
hnite-size scaling relation, valid at the stationary state. A certain quantity Q depends 
on the linear size of the system L and on the deviation e = p — pc according to 

Q(£,L) = (1) 

where is the critical exponents related to the spatial correlation length and x is 
the critical exponent related to Q in the thermodynamic limit, that is, Q ~ when 
L ^ oo. The second type is a time dependent scaling relation, valid for an inhnite 
system, and given by 

Q(£,f) = (2) 

where is the critical exponents related to the time correlation length. When f —)■ cxo, 
we get the same behavior Q ~ At the critical point, relations ([I]) and ([2]) predict 
the following scaling forms 


Q^ty, 


( 3 ) 
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respectively, where y = x/i>\<^. 

The behavior of the SEI model is characterized by a set of quantities. We dehne the 
average Nj = (n/) where nj is the number of infected sites and the average Nsi = (nsi) 
in the number of active pair of sites dehned as the number nsi of pairs of neighboring 
sites of type SI. Another relevant quantity is the surviving probability P{t), dehned as 
the probability that at time t the system is active, that is, as long nsi 7 ^ 0. Starting 
from one infected site in a inhnite lattice full of susceptible, we expect the following 
asymptotic time behavior, at the critical point. 


Nl{t)r~.t\ 

(4) 

Nsi{t) ~ f®. 

(5) 

P{t) ~ t-\ 

(6) 


(7) 


where f is the spatial correlation length and 

If several trials are carried out, some survive up to time t, some do not. In the 
limit f —)■ 00 , the cluster of a surving trial will be identihed with the inhnite percolating 
cluster. Let be the variable that counts the number of infected sites in the surviving 
trials. Denoting by M and Mg the number of trials and the number of surviving trials, 
respectively, and by ngi the number of infected sites in the i-th surviving trial, then the 
average Ng of ng over the surving trials is given by 



The quantity between parentheses, which we denote by Nf is the average number of 
ng over all trials because = 0 for a nonsurviving trial. Taking into account that 
P = Mg/M then Ng = Nj/P. Now, the quantity Nj has a behavior similar to the 
average Ni of ni over all trials. In fact, we found numerically that Nj is proportional to 
Ni for large enough t. Therefore, for large enough t we may write Ng ~ Ni/P. Taking 
into account relations ( 0 ]) and ([ 6 ]), 

Ng{t)r^t^^^. (9) 

The number of infected sites in surviving trials inside a region of linear size f is 
proportional to the order parameter P. This amounts to say that the ratio Ng/^'^ is 
proportional to P. Writing thus ~ P, and taking into account relations o. o 

and (|9]), we reach the following exponent relation 

z{ri + 26) = d. ( 10 ) 

At the critical point a surviving trial, which makes up the percolating cluster, has 
a fractal structure with a fractal dimension dp dehned by 


( 11 ) 
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Figure 1. Static properties of the SEI model from numerical simulations on a square 
lattice, (a) Order parameter P versus p for several values of L indicated. The inset 
is a data collapse showing y = versus x = where e = p — Pc- (b) The 

quantity UP versus p for several values of L indicated. The inset is a data collapse 
showing UP versus x. 


Comparing ([7]) and ([9]) with (HU, we see that the fractal dimension dp oi the critical 
cluster is related to rj and 6 by 

z{rj + 5) = dp- ( 12 ) 


The exponents 9 and rj are connected by the relation 
rj = l + 9, 


(13) 


which can be understood by observing that the rate of increase in the average number 
of infected sites is proportional to the number of active pair of sites, that is, 




(14) 


where k is the number of neighbors. Replacing (|5]) into ffTTj) . it follows that Nj tl+B 
from which we hnd relation ([13]). We should remark that in models belonging to direct 
percolation (DP) universality class [12], W ~ Npi so that the exponents rj and 6 
coincide. In the case of models belonging to dynamical percolation universality class, 
such as the SIR and SEI models, they are distinct and are related by f[T2]) . 

When t ^ oo, that is, in the stationary state, the surviving probability, which is 
identihed with the order parameter, behaves around the critical point as 


P~£^ 


(15) 


and the exponent ft is related to 6 by 


6u\\ = j3, (16) 

which is equivalent to 6z = fi/vp- Taking into account this relation and comparing 
relations (1121) and (ITO]) . we hnd pz = where 7 = dup — 2ft, and (d — dp)h'± = ft. 
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Figure 2. Number of infected sites Ni, number of active pair of sites Nsi, surviving 
probability P and spreading of the infected sites as a function of time t in a double¬ 
log plot. The slopes of the curves for large enough time give, respectively, 77 , 9, —S and 
2/z. 


The density p of infected sites is eqnal to to Nj/N. Bearing in mind that Nj = N^P 
and that Ng = PN, we conclnde that p so that 

p ~ . (17) 

Let ns dehne the qnantity p 2 = {n'j)/N^. Taking into acconnt that (nf) = (nj)/P and 
that {nl)/N‘^ P^ aronnd the critical point, we may conclnde that p 2 P3, or 

P2 ~ (18) 

The qnantity U, dehned as the ratio U = {nj)/(nj)^ = P 2 Ip^ behaves as 

U ~ e-^. (19) 

Scaling relations can be written for the qnantities dehned above. For instance, the 
order parameter obeys the scaling relation 

P(£, L) = (20) 

From the above relations for P and P, it follows that the prodnct UP approaches a 
constant at the critical point and obeys the relation [IT] 

UP = ^2i£‘'^L), ( 21 ) 

so that, at the critical point, e = 0, a plot of UP versns p for several valnes of L will 
cross at the same point, the critical point. 

4. Simulations 

We performed nnmerical simnlation on a sqnare lattice by nsing the algorithms explained 
above. The stationary properties were obtained on lattices of several sizes, and for 
several valnes of the parameter p. The order parameter P was obtained as follows ig. 
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We perform several runs starting from an infected individual placed in the center of a 
finite lattice full of susceptible individuals. The quantity P{p, L) is the fraction of runs 
such that an infected individual reaches the border of the lattice of linear size L. In 
hgure [1^ we show the order parameter P versus p for several values of L. The inset 
of hgure [IK shows a data colapse according to the scaling form fl20l) . To get the data 
colapse we used the exact values of the percolation exponents in two dimensions | 22 j : 
ft = 5/36 and = 4/3. We used also the numerical value of p at the critical point on 
a square lattice, Pc = 0.59274606(5) obtained from Monte Carlo simulations [23] . 

In hgure [T]d, we show the product UP as a function of p. As expected from scaling 
relation fl?Ill the curves for distinct L cross at the critical point. The value of UP at pc 
is estimated to be {UP)c = 1.01658(1), which is a linear extrapolation in 1 /L^ obtained 
from the values 1.016598, 1.016585, 1.016582 of this quantity for L = 256,512,1024 
respectively. The quantity UP at pc is a universal quantity and may be compared to 
the values 1.0167(1) obtained for percolation and SIR models [H]. The inset of hgure 
UK shows the data collapse according to (ETj). Again we use the exact value = 4/3. 

To get the dynamic exponents dehned by relations O, dSD, ©, dB, we performed 
time-dependent Monte Carlo simulations at the critical point pc = 0.59274606(5) [23] . 
To get the exponents rj, 9 and 6, we have estimated the number of infected sites Nj, the 
number of active pair of sites Nsi and the surviving probability P. The exponent z was 
obtained by calculation the spreading of the infected sites [7], dehned by 

“As(22) 

i 

where r* is the distance of site i to the origin and rji takes the value 1 when site i is 
occupied by an infected individual and zero otherwise. At the critical point 

~ (23) 

Figure [2] shows the quantities Nj, Nsi, P and R^ as a function of time. Each curve 
corresponds to an average over 10 ® runs obtained on a square lattice of linear size 
L = 2 ^®. At each run, the lattice was full of susceptible sites except the central site, 
which is in an infected state. Up to the maximum time shown in hgure [2]no infected site 
reached the border of the lattice, which amounts to say that the results shown in hgure 
HI are valid for an inhnite lattice. We remark that the quantity Nj was also calculated 
and we found cNj = Nj for large enough t, as we have assume above, and that at the 
critical point c = 1.025(5). 

To estimate the exponent y for a given quantity Q ^ C, we used a correction to 
scaling of the form [7] 

Q{t) = U{ci + C2t~^). (24) 

The htting of this form to the data points gives the exponent y and y in addition to 
the constants ci and C 2 . The htting was done within a certain interval of time starting 
at time R and ending at time R. We have found that the best httings give y around 1, 
but the actual values found for the exponent y are not too sensitive to the value of the 
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Table 1. Critical exponents obtained by fitting the scaling form (|4l) to the curves in 
figure [2] within the time interval between to and ti. The last row gives the average of 
the values in the previous rows, together with the errors in the last digit. 


Into 

Inti 

p 

9 

5 


4.5 

9.5 

1.58462 

0.58444 

0.09217 

1.13096 

5.5 

9.5 

1.58457 

0.58438 

0.09208 

1.13087 

6.5 

9.5 

1.58454 

0.58442 

0.09196 

1.13083 



1.5846(2) 

0.5844(2) 

0.0921(3) 

1.1309(3) 


exponent /t. We used three values of to and the same value of ti. As seen in table [U, the 
three values of a given exponent are distinct but very similar allowing the estimation of 
the statistical error shown in the last row together with the average. 

The estimated exponents can be checked by using relations pz = and 

6z = From the values of p and we get pz = 1.7920(7) which should be 

compared with the exact result = 43/24 = 1.791666... [22]. From the values 

of 6 and 2 ; we get 6z = 0.1042(4) which should be compared with the exact result 
= 5/48 = 0.104166... [22|. We remark that our results are in agreement with 
results for the exponents of models belonging to the universality class of dynamic 
percolation [ISl EO]. Notice that the critical exponents p and 6 satisfy, within the 
statistical error, the relation ffT3l) . p = 1 -\- 9. 

5. Discussion 

We have determined the critical properties of the SEI model on a square lattice by 
numerical simulations. The model is an example of model belonging to the universality 
class of dynamical percolation. The stationary properties are shown to be exactly 
the same as those of isotropic percolation, so that the static exponents are the same 
as percolation, which in two dimensions are known exactly. This is not the case 
of the dynamic exponents which we have calculated here by using a time-dependent 
simulations. Our numerical estimation, within the statistical errors, are in good 
agreement with the known relations among the exponents. An important feature that 
distinguish this model from models belonging to the DP universality class na, such as 
the susceptible-infected-removed-susceptible (SIRS) [T3|, rests on the relation between 
the number of infected sites and the number of active pairs of sites. In the SIRS model 
these two quantities are closed related leading to the identihcation of the exponents p 
and 9. In the SEI model they are distinct leading to the relation p = 1 -\- 9, which was 
conhrmed by our numerical simulations. 
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